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For a special class of bipartite states we calculate explicitly the asymptotic relative entropy of 
entanglement E^ with respect to states having a positive partial transpose (PPT) . This quantity is 
an upper bound to distillable entanglement. The states considered are invariant under rotations of 
the form 0®0, where O is any orthogonal matrix. We show that in this case E^ is equal to another 
upper bound on distillable entanglement, constructed by Rains. To perform these calculations, we 
have introduced a number of new results that are interesting in their own right: (i) the Rains bound 
is convex and continuous; (ii) under some weak assumption, the Rains bound is an upper bound to 
E^; (iii) for states for which the relative entropy of entanglement En is additive, the Rains bound 
is equal to Er. 

PACS numbers: 03.67.-a, 03.67.Hk 



I. INTRODUCTION 

In spite of the impressive recent progress in the the- 
ory of entanglement [1], many fundamental questions or 
challenges still remain open. One of these issues is to de- 
cide whether a given state is entangled or not. Another 
question is to find criteria for the distillability of a state, 
i.e. whether pure state entanglement can be recovered 
from the original state by means of local operations and 
classical information exchange. 

Since entangled states are a resource in many basic 
protocols in quantum computation and quantum commu- 
nication, a need has emerged to quantify entanglement. 
This leads to more advanced challenges: how much en- 
tanglement is needed to create a given state and how 
much entanglement can be recovered? 

Since these questions lead to very high dimensional op- 
timization problems, it is often helpful or even inevitable 
to restrict oneself to states exhibiting a very high sym- 
metry. The two most common one-parameter families of 
symmetric states are the so-called 'Werner' states [2] and 
the 'Isotropic' states, which are related to one another via 
the partial transposition operation. A larger set of sym- 
metric states, containing these two sets as special cases, 
are the OO-invariant states, which are the states consid- 
ered in this paper. 

So far it is not known how to calculate distillation rates 
for arbitrary states, and even for symmetric states this 
optimization seems to be intractable. One possible way 
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to partially circumvent this problem is to calculate good 
bounds for the distillation rates. A well-known upper 
bound for the distillable entanglement is the relative en- 
tropy of entanglement [3], which is itself defined as an 
optimization: 

E R {p) = inf S{p\\a). 

In this formula, S(p\\a) — Tr(plogp — ploga) is the 
relative entropy (the quantum mechanical analog of the 
Kullback-Lcibler divergence) and the minimum is taken 
over all states a in the convex set T>. The relative entropy 
between two states is a measure of distinguishability and 
can intuitively be regarded as a kind of distance measure, 
although it violates most of the axioms that are required 
of a distance measure [3] . In the originally proposed def- 
inition of the relative entropy of entanglement, T> is the 
set of separable states, so that the E R {p) expresses the 
minimal distinguishability between the given state and 
all possible separable states. When using the Er as an 
upper bound to distillability, however, it is fruitful to en- 
large the set D to the set of states with positive partial 
transpose (the PPT states) [4] . The corresponding mini- 
mal relative entropy, the relative entropy of entanglement 
with respect to PPT states (REEP), is generally smaller 
than the (separability) relative entropy of entanglement 
while it still is an upper bound to distillability; this is 
so because all PPT states have distillability zero. Hence, 
the REEP is a sharper bound on the distillability than 
the separability relent. This enlargement of V has the ad- 
ditional benefit that the set of PPT states is much easier 
to characterize than the set of separable states, for which 
no general operational membership criterion exists. 

Nevertheless, neither for the REEP nor for the relative 
entropy of entanglement is there a general solution known 
of the optimization problem for arbitrary states, not even 
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for the otherwise simple case of two qubits. However, the 
calculations become tractable when restricting oneself to 
symmetric states. 

Contrary to earlier conjectures, neither REEP nor the 
relative entropy of entanglement is additive, i.e. Er(pi ® 
p-i) < Er(pi) + Eji{p 2 ) is a strict inequality for some 
states. It is expected, however, that this non-additivity 
becomes less severe for the asymptotic relative entropy of 
entanglement with respect to PPT states(AREF,P) , which 
is defined as the regularisation 

E%{p) = lim -E R (p^), 

n—>oc n 

and which at the same time provides yet a sharper bound 
to distillable entanglement. 

The calculation of the AREEP has first been done on 
Werner states [10], showing that the asymptotic value can 
be a good deal smaller than the single-copy value. Sur- 
prisingly, it turns out that on Werner states the AREEP 
is equal to another upper bound on distillability, the so- 
called Rains bound [5] 

R(p) =infS(H|a)+logTr|a T2 |. (1) 

(7 

One of the things we will show in this paper is that 
this equality remains valid over the larger class of 00- 
invariant states. 

To calculate the AREEP on OO-invariant states in a 
relatively simple way, we will make use of four ingredi- 
ents: 

• First of all, REEP is additive on a large part of 
the state space. This will be discussed in Sec. 
II. For this additive region, the calculation of the 
AREEP is trivial, as the (single-copy) REEP for 
OO-invariant states has been calculated before. 

• We will make use of the convexity of the AREEP 
(recollected in Sec. Ill) and of the Rains bound 
(proven in Sec. IV). In Sec. Ill we use this convexity 
to define the "minimal convex extension" of the 
AREEP from the additive areas to the full state 
space. 

• In Sec. V we will present a close connection between 
the Rains Bound R(p) and the AREEP. We will 
establish an upper bound to the AREEP that will 
turn out to be tight on OO-invariant states. 

• In Sec. VI B we will recall the basic properties of 
OO-invariant states resulting from their symmetry. 
It is exactly this symmetry that makes the calcula- 
tion feasible. 

Using these results, we will give a complete calculation of 
the AREEP of OO-invariant states in Sec. VI and prove 
that this quantity is equal to the Rains bound for these 
states. We will summarise the results of the paper in Sec. 
VII and state a number of open problems. 



II. ADDITIVITY OF RELATIVE ENTROPY OF 
ENTANGLEMENT 

The additivity of the REEP was a folk conjecture, sup- 
ported by various numerical calculations and analytical 
case studies. Nevertheless, it turned out to be wrong 
[6]. The misleading numerical result can be explained 
in hindsight by the fact that, indeed, in great parts of 
the state space the REEP is perfectly additive; the non- 
additive regions seem to be negligible in size compared 
to the whole state space. 

The following Lemma of Rains [4] can be utilized to 
pinpoint regions where the REEP is additive. 

Lemma 1 (Rains- Additivity) Let p be a state and a 
a PPT state, such that Er(p) = S(p\\a) and [p,a] = 0. 
// the condition 

|(pO T2 l<* (2) 

holds, then the REEP is weakly additive on p, i.e., 
Eft(p) = Er{p). If it satisfies the stronger condition 

< (p^f 2 < t (3) 

then REEP is strongly additive, i.e., Er{p®t) — Er(p) + 
Er(t) holds for an arbitrary state r. 

Knowing the optimal a for a given state p, it is straight- 
forward to check condition (2). Checking the additivity 
therefore only requires one to calculate the REEP. 

III. CONVEXITY OF THE ASYMPTOTIC 
RELENT 

By definition, the asymptotic version of a given quan- 
tity inherits most of the important properties directly 
from its single-copy "parent" quantity. One such prop- 
erty, which will turn out to be very helpful to calculate 
the AREEP, is convexity. The REEP itself is known 
to be convex, but it is not obvious that quantities of the 
form E n (p) :— E(p® n )/n should be convex functions in p 
too and, in fact, this does not hold in general. Although 
convexity might not hold for finite n, for the REEP it 
becomes valid again in the asymptotic limit. 

Lemma 2 [9] Let E be a positive, subadditive, convex 
and tens or- commutative functional on the density ma- 
trices of a Hilbert space. Then the asymptotic measure 
E°°(p) :— linin^oo ^E(p® n ) exists and is convex and 
subadditive. 

In the first calculation of the AREEP [10] great effort 
was necessary to construct a lower bound to AREEP. 
Utilizing the convexity we arc now able to do this in a 
much simpler way. Indeed, for any convex (differentiable) 
function /, a lower bound to / is given by any of its 
tangent planes 

f(x)>f(y) + Vf(y)(x-y). 
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Given an open subset D where the function / is known, 
we can define the "minimal convex extension" of the func- 
tion by 

f(x) = sup f(y) + V/(y)(x - y). 

Note that / is equal to / on D. Furthermore, / is smaller 
or equal than any convex function that equals / on D. 
As a maximum over affine functions it is itself convex. 

To make this bound a good candidate for an estima- 
tion to the AREEP, we need to know the AREEP on 
a sufficiently large part of the state space. In fact the 
AREEP is easy to calculate on PPT states, where it is 
simply zero. But this is obviously too trivial a result, 
because this gives a lower bound equal to zero on the 
whole state space. The next greater set for which we can 
easily calculate the AREEP is the set of states where 
is additive. A subset of this set can be found using the 
Lemma of Rains. It will turn out that this subset is large 
enough to yield a bound that equals (at least for 
OO-invariant states). 

IV. CONVEXITY AND CONTINUITY OF THE 
RAINS BOUND 

Although the function that is to be minimized in Rains' 
bound, S , (p||cr)+logTr \a T2 1, is not convex in a over state 
space, the minimum itself turns out to be convex in p. We 
prove this by first showing that the minimization problem 
in the calculation of the Rains bound can be converted 
to a convex problem. 

To begin with, we can add a third term to the function 
to be minimized, namely — logTr[cr], because this term is 
zero anyway. Secondly, we can enlarge the set over which 
one has to minimize from the set of normalized states to 
the set S = {s > 0,Tr[s] < 1}. This is so because the 
sum of the first two terms is independent of Tr[cr] and 
the third one monotonously decreases with Tr[c]; hence, 
the minimal value must be found on the boundary of S 
corresponding with Tr[cr] = 1 and is, therefore, equal to 
the original minimum. The second and third term can 
now be absorbed in the first term: S(p\ \a) +log Tr \a T2 \ — 
logTrcr = 5(p||cr(Trcr/Tr|cr T2 |)). Defining 

T = ( r(Tr ( T/Tr|cr T2 |), 

it is easy to check that a G S if and only if r G T = 
{t > 0, Tr \t T2 1 < 1}. Hence, the calculation of the Rains 
bound has been transformed to the minimization prob- 
lem 

R(p) = minS(p\\r). 

The importance of this transformation stems from the 
fact that the resulting optimization problem is a so-called 
convex optimization problem: the function to be mini- 
mized is now convex in t, while the set over which the 



minimization is performed is still convex. The latter 
statement follows directly from the convexity of the nega- 
tivity. Indeed, if t\ and T2 are in T, then they are positive 
and have negativity < 1. Hence, any convex combination 
of Ti and T2 is positive and has negativity < 1 as well, 
and, therefore, belongs to the set T. 

It is now easy to prove continuity and convexity of the 
Rains bound itself. Continuity follows by noting that the 
proof of continuity of the quantity inf CTe x> S^pH^) in [11], 
where T> is a compact convex set of normalized states 
containing the maximally mixed state, does actually not 
depend on the trace of the various a in T>. Hence, the 
theorem is also true for convex sets T> containing non- 
normalized states, and, specifically, for the set T. 

Convexity is also proven in the standard way, as has 
been done for En [3] . The standard proof again depends 
only on the convexity of the feasible set and not on the 
normalization of the states it contains. 

In this way we have proven: 

Lemma 3 The calculation of the Rains bound can be re- 
formulated as a convex minimization problem: 

R{p) = min{5'(p||T) : r > 0,Tr|r T2 | < 1}. 

The Rains bound itself is a continuous and convex func- 
tion of p. 

V. RELATION BETWEEN RAINS' BOUND 
AND THE AREEP 

The results of the calculation of the AREEP on Werner 
states suggests [7] that this quantity might be connected 
with the quantity (1) defined by Rains, and, moreover, 
that there are connections between the minimizing a in 
Rains' formula and the asymptotic PPT state a appear- 
ing in Ej£ . Indeed, it turns out that one can give a simple 
relation between these two quantities, if we require as an 
additional restriction that a in (1) satisfies \a T2 \ T2 > 0. 
If the restriction does not hold the Lemma might still be 
true, but we have not been able to prove this. 

Lemma 4 An upper bound for the AREEP is given by 

R'(p) := infS(p||a)+log(Tr|<7 T2 |) 

> E%(p), (4) 

where the asterisk means that the infimum is to be taken 
over all states a satisfying 

\<r T2 \ T2 >0. (5) 

We will refer to the quantity R'(p) as the modified Rains 
bound. The set of states satisfying condition (5) is easily 
seen to be convex, so that the modified Rains bound is 
also continuous and convex. 

Proof: It can easily be seen that the Lemma is valid if we 
restrict a to be a PPT state, since then the second term 
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in (4) vanishes and we get the trivial inequality E^ (p) < 
E R (p). This means that we can restrict ourselves to the 
case where a is a non- PPT state, i.e. Tr|cr T2 | > 1. 

Let a be an arbitrary non-PPT state such that a :— 
\ a T2\T 2 > o, then 



a n = 



(Ira) 



is a PPT-state. Taking this PPT state as a trial state in 
the optimization for the AREEP, we get 



E R (p® n ) < S{p® n \\a n ) = S p® n 



_|_ Q-® n 

1 + (Tr u) n 



^ S \P "l + (Tr^ 
= n5(p||a)+log(l + (Tra) 



(6) 



In (6) we have used the fact that the relative entropy is 
operator anti-monotone in its second argument (Corol- 
lary 5.12 of [8]), i.e. S{p\\a + r) < S(p\\a) for positive r. 
Taking the limit n — > oo and using Tr a > 1 we get 



hm -E R (p® n ) 

n— »oo XI 



< lim S(p\\a) + 



log(l + (Tr ( r) 11 ) 



S(p\\a) +logTra. 



(7) 



In order to get the best bound, we take the minimum 
over all feasible states a in equation (7), giving 



££(p)<infS(p||a)+logTr|a 



T 2 



where the infimum is taken over all states a satisfying 
|ct T2 | T2 >0. □ 



It is easy to see that for PPT states a, \a T2 \ T2 > 0. 
Hence, the feasible set in the minimization of E R is a 
subset of the one for R', which is again a subset of the 
one for R. Therefore, we have the inequalities 

R(p) < R\p) < E R (p). 

We also have the following Theorem: 

Theorem 5 For E R -additive states p (i.e. E R (p) = 
E R c (p)), the Rains bound is equal to the AREEP and 
is additive. 

Proof: We have, in general, R'(p) < Eji(p). On the other 
hand, for additive states E R (p) = E R °(p), and E^(p) < 
R'(p) by the Lemma 4. Therefore, R {p) = E R (p) = 
E^ (p) for all additive p. This also implies that the PPT 
state a that is optimal for E R is also optimal for R' . 

To show that R is also equal to E R , we need to show 
that this a is optimal for R as well. We use the refor- 
mulation of the Rains bound as a convex minimization 



problem R(p) = mh\ T {S{p\ |r) : Tr|r T2 | < 1}. For the 
modified Rains bound, we have the additional restriction 
on the feasible set that |t T2 | T2 > 0. For clarity, let us 
write r for the optimal r for R and r' for the optimal 
one for R'. We have to show that r = r', i.e. that r is in 
the set for which |t T2 | T2 > 0. 

Suppose r were outside this set, then, following a gen- 
eral property of convex optimization problems, r' would 
have to be on the boundary of the set, i.e. \t' T2 \ T2 
would have to be positive and rank-deficient. On the 
other hand, we already showed that the optimal a' for 
R' for additive p must be PPT, so that r' = a' and 
\a' T2 \ T2 = a' . Therefore, the rank-deficiency of \t' T2 \ T2 
implies that a' itself should be rank-deficient. However, 
if p is not itself rank-deficient, then this cannot be, be- 
cause a' appears as second argument in the relative en- 
tropy and would then give an infinite relative entropy, 
contrary to the statement that a' actually minimises it. 
This proves that R'(p) = R(p) for full-rank, additive p. 
By continuity of the Rains bound this must then also 
hold for rank-deficient p. 

Additivity of R for _E^-additive states follows by reg- 
ularising both sides of the equality R(p) — E R °(p), and 
noting that the right-hand side does not change. □ 



We have introduced the operation a i— > \a T2 \ T2 as a 
mathematical tool, and we doubt whether it has any real 
physical significance (as was the case for the partial trans- 
pose). Nevertheless, its usefulness is apparent from the 
above Lemma. A natural question to ask is whether there 
really are states a for which \(j T2 \ T2 is not positive. We 
call states like this binegative states. If they would not ex- 
ist, then the modified Rains bound would just be equal to 
the original Rains bound. We have performed numerical 
investigations that have shown that, indeed, binegative 
states exist, provided the dimensions of the system are 
higher than 2x2. For 2x2 systems, extensive calcula- 
tions failed to produce binegative states, which suggests 
they might not exist in such systems. For higher dimen- 
sions, binegative states have been produced, and they 
always appear to be located close to the boundary of 
state space, i.e. have a smallest eigenvalue which is very 
small. In the present setting, this is good news, because 
it implies that the modified Rains bound will typically 
be close to the original Rains bound. 

As one of the few exact results on the existence of 
binegative states, we have been able to prove that pure 
states are never binegative: 



Lemma 6 For any pure state ip, 



T 2 



T 2 



> 0. 



Proof: Let ip have a Schmidt decomposition \p = 
YsiKuit&Vi, then I^XV'F 2 = ^i^j\ u i)( u j\ ® \ v i)( v j\ T 
and, exploiting the orthogonality of the vectors m and of 
the vectors Vj, 



( ^2 ^j^kh\ui)(uj\\u k )(ui\® 

i,j,k,l 
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(k)Hk>H) T ) 1/2 

= (x;(^) 2 i«iX«ii®i«iX^-r) 1/2 . 

since only the terms with i = I and j = k survive. Again 
by orthogonality, taking the square root amounts to re- 
moving the square on the factor (XiXj) 2 . Now, one clearly 
sees that the resulting expression corresponds to a prod- 
uct state; hence, the partial transpose is still a state, 
which proves that ip is not binegative. □ 



In Section VI B we will show that no OO-invariant 
state is binegative either. We will see that condition (5) 
will be fulfilled for the states we are considering in this 
paper. Therefore, we will henceforth make no distinction 
anymore between R and R' . 



VI. OO-INVARIANT STATES 

We will now apply the tools obtained in the previous 
sections to the complete calculation of the AREEP of 
OO-invariant states. 



In the last step we calculate the tangent to the REEP 
in the point p = | + i , which gives us the minimal convex 
extension for all states with p > \ + ^ . It turns out that 
this minimal extension touches the Rains bound again 
in the point p = 1. This is sufficient to prove that the 
minimal convex extension is equal to E^ everywhere. 
Indeed, by the convexity of E^ the tangent yields a lower 
bound and, furthermore, also implies that the tangent is 
an upper bound between p = 1 and p = i + i , because 
at the end-points it equals E^. 

In fact, for Werner states, the same result can eas- 
ily be obtained by the observation that the Rains bound 
and the minimal convex extension are equal on the whole 
range of p. But for OO-invariant states the task to prove 
equality of these two quantities will become quite dif- 
ficult. Fortunately, we can restrict ourselves to prove 
equality only on the border of the state space as this will 
be sufficient for the calculation. Equality of the Rains 
bound and on the whole state space will follow au- 
tomatically from the convexity of both quantities. 

We will now turn to the calculation for the OO- 
invariant states. 

B. Using symmetries 



A. Calculating the AREP on Werner states 

To illustrate how the calculation of the AREEP on 
OO-invariant states will proceed, we apply the method 
first on Werner states, reproducing the results of [10]. 

Werner states can be written as 

P P 
p(p)=p— + (l-p) — , 
r_ r + 

where denotes the normalized projection onto 

the symmetric (antisymmetric) subspace of dimension 

j 2 I j 

r± = and p is a real parameter ranging from 

to 1. 

First of all, we need to know E R on these states. All 
states with P <\ are PPT and, therefore, have both Er 
and E^ equal to zero. For all non-PPT Werner states 
p > \ , the minimizing PPT-state is the state with p =\- 
Knowing this state, we can easily write down the REEP 
for all Werner states. To calculate the AREEP we use 
the three steps introduced in the previous three Sections. 

In the first step we use the lemma of Rains and check 
the additivity condition (2). An easy and straightforward 
calculation leads to the result that all Werner states sat- 
isfying p < i + 2 arc additive and, therefore, have Er 
equal to E^. 

In the second step we calculate the Rains bound for 
Werner states. Due to the high symmetry this is an easy 
task, already done by Rains [5]. In fact, we do not need to 
compute the Rains bound for all states. For our purposes, 
we will only need the Rains bound for p = 1. 



The class of states we want to look at commute with 
all unitaries of the form O ® O, where O is an orthogo- 
nal matrix. These so-called OO-invariant states lie in the 
commutant G' of the group G — {O ® O}. The commu- 
tant is spanned by three operators, the identity operator 
1, the Flip operator F defined as the unique operator for 
which Wip<3cf) = (j><S)tp for all vectors tp and <fi, and the un- 
normalized projection on the maximally entangled state 

^ = J2ij N)0'il — ^I^X^Ii nere 7 d is the dimension of 
either subsystem. Every operator contained in this com- 
mutant can be written as a linear combination of these 
three operators. To be a proper state such an operator 
has to fulfill the two additional constraints of positivity 
and normalization. 

As coordinates parameterizing the OO-invariant 
states, we choose the expectation values of the three op- 
erators 1,F and F in the given state. The expectation 
value of the identity, (l) p , gives us just the normaliza- 
tion, so we are left with the two free parameters / := (F) 

and / := (F)^. For future reference, we collect the basic 
formulae here for performing calculations in this repre- 
sentation. 

The traces of the basis operators are given by 

Tr[l] = d 2 
Tr[F] = d 
Tr[F] = d. 

The inner products between them are easily calculated 
from the relations 

F 2 = 1 
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FF = FF = F 
F 2 = d¥. 

From this basis {1,F, F}, an orthogonal basis of projec- 
tors can be constructed. The operator F is not posi- 
tive and can be written as F = F + — F_; here F + and 
F_ denote the positive and negative part of F, respec- 
tively, and are defined by the equations X = X + — X- , 
\X\ = X + + X- (note that both the positive and nega- 
tive part are positive by this definition). Since F 2 = 1, 
F + + F_ = 1, and F_ = (1 - F)/2. Furthermore, as 
FF = F, F < F+. Therefore, the following operators 
form an orthogonal set of projectors and add up to the 
identity: 

U = F/d 

V = (l-F)/2 

W = (l+F)/2-F/d. 

The traces of these projectors are 

Tr[U] = 1 
Tr[V] = d(d -l)/2 
Tr[W] = (d + 2)(d- l)/2. 

The original basis is related to the orthogonal one by 

1 = U + V + W 
F = U-V + W 
F = dU. 

For a general OO-invariant p, we write 

p = at + b¥ + cF. 

The relation between the coefficients a, b, c and / and / 
is given by 



1 — 1 




d 


1 


1 


/ 


= d 


1 


d 


1 


/. 




1 


1 


d 



and, inversely, by 



1 



1 



-1 -1 




1 


d + 1 -1 




/ 


-1 d + 1 




/ 



d(d- l)(d + 2) 
In terms of the orthonormal basis, p can be written as 
1-/ t r d + df-2f 



?= f d u 



-V 



d(d- 1) d(d- l)(d + 2) 
Positivity of p thus amounts to the conditions 

< / 
/ < 1 

/ < d(l + /)/2. 



W. (8) 



The representation of the partial transpose of p is very 
easy, since F and F are just each other's partial trans- 
pose. Hence, the partial transpose of p is obtained by 
swapping F and F. In the basis { 11 . 7. 7} . taking the 
partial transpose corresponds, therefore, to interchang- 
ing the parameters / and /. The partial transposes of 
the projectors U, V and W are easily calculated to be 



U l2 
W T2 



-(U-V + W) 



= ( 



1 + d 1, A 
— -d^ + (2 + 



From these formulae one can see that the set of OO- 
invariant states constitutes a triangle in the (/, /) pa- 
rameter space, as plotted in Figure 1. Taking the partial 
transpose amounts to taking the mirror image around the 
line / = /. Therefore, the set of PPT states are those 
contained in the grey square < /, / < 1 in Figure 1 . 

What will make the calculation of the REEP easy for 
these OO-invariant states is the existence of a 'twirl' op- 
eration [2], a projection operation T that maps an ar- 
bitrary state p to an OO-invariant state T(p) and that 
preserves PPT-ness, i.e., that maps every PPT state to 
an OO-invariant PPT state. Since 

S(p\\a)>S(T(p)\\T(a)) 

this guarantees that the minimum relative entropy for an 
OO-invariant state is attained on another OO-invariant 
PPT state [4, 6]. Hence, we can reduce the very high- 
dimensional optimization problem to an optimization in 
our two-dimensional OO-invariant state space. This op- 
timization has been done [6] and the minimizing PPT 
states are as follows. Let a state p be determined by the 
expectation values (F)^ = / and (F)^ = /. Similarly, let 
the expectation values in the optimizing PPT state a be 
given by (F) = s and (F) = s. Then the following ta- 
ble gives the expressions for s and s, depending on which 
region the state p is in: 



Region 


s 


s 


A 


i+(d-i)f-f 


1 




d-f 




B 





/ 

1+/ 


C 





1 



To end this section, we give the formulas for the rel- 
ative entropy and the negativity of OO-invariant states. 
Let the states p and a be determined by their expec- 
tation values /, / and s, s, respectively. Using the state 
representation (8), in the orthogonal basis {U, V, W}, the 
relative entropy of p w.r.t. a is given by 



S(p\\a) = ^log(£)Tr£/ 



W 
d(d-l) 



log 



£3 



TrV 
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FIG. 1: State space of OO-invariant states (case d = 3). 
These states are parameterized by the two parameters / = 
(F) and / = (F). The outer triangle represents the values 
corresponding to states (positivity). The grey area is the set 
of PPT OO-invariant states. The region of non-PPT states is 
subdivided further in the three triangular regions labeled A, 
B and C. For each of these regions the optimal a appearing 
in the definition of the REEP is of a different form. 



d(d-l)(d + 2) 
l-f 



f, f 
-jlog- 

d s 



2 

df - 2/ 



log 



d + ds- 2s 
l-f 



2d 



log 



1 - s 
d + df-2f 
d + ds — 2s 



(9) 



Recollecting that taking the partial transpose corre- 
sponds to interchanging s and s, the negativity of a is 
given by 



Trl^ 2 ! = 



TrU 



1 



d(d-l) 



TtV 



+ 


d + 


ds - 2s 


Tr W 




d(d- 


-l)(d + 2) 




+ ' 1 


-s\ \d- 


-ds- 2s\ 


d 




2 


2d 



(10) 



The positivity condition on a implies that the absolute 
value sign on the third term is superfluous. 



In a similar way, we can show that for any OO- 
invariant state a, the operator \a T2 \ T2 is a state again, as 
we had promised. Indeed, 



,T a |T 2 



U 



T-2 



d(d-l 
d + ds — 2s 



d(d-l)(d + 2) 



V 



W 1 



T 2 



An easy but somewhat lengthy calculation shows that 
this expression can be rewritten in terms of U , V and W 
with positive coefficients. 



C. Additive Areas 

In the first step we want to identify the areas within 
the state space triangle where the REEP is additive. 

Lemma 7 Er(p) is additive for all OO-states satisfying 
(F) > =f and (F) < 3 - $ + (d - !)( F >- 

Proof: Utilizing Lemma 1 we only have to check condition 
(2) for every OO-invariant state p and the corresponding 
optimal PPT-states a. 
directly given by 



In the {U, V, VFj-basis, pa is 



with 



pa 1 = uU + vV + wW, 



s 

l-f 
1 - s 

d + df-2f 



w 



d + ds — 2s 



In order to perform the partial transpose, we replace 
U, V, W by their partial transposes and express them in 
the original U, V, W again. This yields 



(pa' 1 ) 12 = {a + bd + c)U + (a- c)V + (o + c)W, 



with 



w + v 



Condition (2) is then satisfied if and only if \a + c + bd\, 
\a + c| and \a — c\ are all < 1. For s and s we have to 
insert the values of the optimal PPT state a, obtained at 
the end of the previous section. 

After a tedious calculation, we get 6 conditions an ad- 
ditive state has to satisfy for each of the three regions A, 
B and C of Figure 1. Fortunately, only two of this total 
of 18 conditions can be violated by expectation values be- 
longing to normalized positive states. In the A region all 
states are additive, in region B we must have / > — 2/d, 
and in region C the condition is / < 3 — A/d + (d — 1)/. 
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These conditions give us the border between the additive 
and lion additive areas. □ 



The additive area for OO-states is plotted in dark grey 
in Figure 2 for the dimension d = 3. States in the light 
grey area fulfill the condition of strong additivity. 

For later use, we have marked some points in the state 
space that will become important in the further calcula- 
tion of the AREEP. The two additivity conditions of the 
Lemma correspond to the boundary line segments CD 
and BC, respectively. 



(F) 



Point 


<F) 


<F) 


A 


-1 





B 


d-4 
d 


d-2 


C 


-2 
d 


d-2 
d 


D 


-2 
d 





E 





1 


X 
Y 


4-Qd+d 2 


d 2 (d-2) 


d(d+2)-4 
-d 2 


d(d+2)-4 
d(d-2) 


d(d+2)-4 


d(d+2)-A 



D. Rains upper bound 



In the second step we want to calculate the Rains 
bound (4) on the OO-invariant state space. All 00- 
invariant states satisfy \a T2 \ T2 > 0, and we therefore re- 
strict the optimization to OO-invariant states a. Since 
we want to use the Rains bound as an upper bound, we 
need not to know that our so restricted a is really the 
optimal one. But due to the high symmetry of the OO- 
states it can easily be shown that the optimum over all 
possible states a is attained on OO-states anyway. 
For additive states we have noted already that Er (p) = 
R(p) so that calculating the REEP directly 
gives the Rains bound. To calculate the Rains bound 
in the non-additive region ABCD, we have to perform 
the minimization explicitly. Let the states p and a be 
determined by their expectation values /, / and s, s, re- 
spectively. Using the formula for the relative entropy of 
p (9) w.r.t. the optimal a for the REEP (see the Table in 
Section VI B) yields the Rains Bound for additive states. 

For non-additive states we have to include the nega- 
tivity of er, given by (10): 



Er°{p) 



Ti\a A 



d 



d + ds -2s 
2d ' 



As we will only use the above formulae for p in the non- 
additive region ABCD, it is immediately clear from Fig- 
ure 2 that the optimal a will have negative s. We can, 
therefore, simplify the formula for the negativity to 



Trk J 



1 



d 2 
max(l, s) - 



d + ds - 2s 
2d 



2.s 
1' 




<F) 



FIG. 2: Additive areas for OO-invariant states (case d = 
3). The state space has been subdivided in three regions. 
According to Rains' Lemma, the states in the light-grey region 
are strongly additive and those in the dark-grey region are 
weakly additive. The region of additivity is delineated by the 
line segments BC and CD. The points A, B, C, D and E are 
defined in the text. 



Because of the 'max' function appearing in this formula, 
we have to consider two cases for a and, in the end, choose 
the solution that gives the smallest value for the Rains 
bound. 

Consider first the case s > 1; then the negativity equals 
s — 2s I d and we have to minimize 



log- 



ics - 2s 



1 

2d 



1 



+(d + d/-2/)lo| 



2/log( + (d-df)logl— 
s s — 1 



d + df-2f 



d + ds — 2s 



over s and s. This function has a single stationary point 
given by 

d 2 -df-2 



(d 2 - 2)/ - df 

-2/ 
(d 2 -2)f-df 



However, the minimum we are looking for is a constrained 
one: the parameters s and s must be expectation values 
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of positive er. On inspection, the positivity conditions 
are never satisfied in the stationary point for any choice 
of /, / corresponding to a positive p. Therefore, the sta- 
tionary point is outside the feasible set (the state trian- 
gle) and the constrained minimum will be found on the 
boundary of the feasible set. This mere fact already rules 
out the present case s > 1 , because we know that the op- 
timal a must be closer to the set of PPT states than p 
itself, in the sense that a should have lower negativity 
than p. Indeed, setting a = p (which is certainly not op- 
timal) in the Rains bound yields a lower value than one 
would get for any a with a larger negativity than p. 

We can, therefore, restrict ourselves to the case s < 1. 
As the negativity is then 1 — 2s/ 'd, the function to be 
minimized is 



R(P) 







B d=5~ 


A 






a\ 




/x 


\ \ Y 
a\ \ 




B d=4 














Y \^ 


B d=3 

X 




-0.8 


-0.6 -0.4 -0.2 


0.2 



log 



d-2i 



1 

2d 



2/ log 4 + (d - df) log 1—^ 

s s — 1 



+(d+d/-2/)lo| 
The stationary point is 



d + df-2f 



d + ds — 2s 



2 + df 
d+2f 

(2 + d)f 
d + 2f ■ 



(11) 

(12) 
(13) 



Again, s and s must be expectation values of positive 
a and we must have that s < 1. It turns out that the 
positivity conditions arc always fulfilled. The condition 
s < 1, on the other hand, is only satisfied for states p on 
or below the line going through points C and Y. There- 
fore, the stationary point is the constrained minimum 
only for states p in the quadrangle AYCD. This leads to 
the solution for AYCD: 



Ra 



YCD 



(P) 



1 



J(l + /)log(d-2)-21ogd 
-(/-l)bg(d + 2)), 



(14) 



which now only depends on the Flip expectation value / 
and is an affine function of /. 

For states p in the remaining triangle CYB, the sta- 
tionary point is outside the feasible set, so that the con- 
strained minimum will lie on the line 1=1. Minimization 
of (11) over s, while fixing s = 1, yields a quite cumber- 
some looking formula. For later use, however, we will 
only need to know the resulting Rains bound on the line 
segment YB. The solution consists of two cases, corre- 
sponding to either solution of a quadratic equation. The 
end result is that, for the states on the segment YX, the 
Rains bound is given by 



Ryx(p) 



1 + / 
2 
1 



logd(l + /) 
/ d(l-f) 



, d{d + 2)-4 , n 
+ log- — ^ log 2. 



(15) 



FIG. 3: Rains bound on the line segment AB (see Figure 
2) in terms of the parameter /, for three different values of 
d = 3, 4, 5. The bound consists here of a linear part (segment 
AY, eq. 14), a curvilinear part (segment YX, eq. 15) and again 
a linear part (segment XB, eq. 16). 



For the states on the segment XB the bound is given by 



Rxb(p) 



1 + / 



Iog(d - 2) 



f 



-1 d 
log — . 

2 S 4 



(16) 



Figure 3 shows the Rains bound along the line segment 
AB, for several dimensions d = 3, 4, 5. 



E. Minimal convex extension 

In this third and final step we calculate the minimal 
convex extension of the additive area. This will turn 
out to be more complicated than in the Werner states 
example. We will look at straight lines, each connecting 
one point on the additivity border with one, well-chosen 
point on the line segment AB. 

The simplest case is the part of the additivity border 
consisting of the line segment CD, because this line lies 
completely in the 'Werner' region, region B in Figure 1, 
where, according to (14), the REEP only depends on the 
Flip expectation value /. So, here, the two-dimensional 
problem is reduced to a one-dimensional one. The REEP 
in the Werner triangle is given by 

£B(/) . log2+ (l±fl log i±/ + <i_fl log i_/. 

As lower bound for the AREEP we get 

> E R (-2/d) + (f + 2/d)^^-\ f= _ 2/d 



/)log 



d-2 
d+2 



Of 

log^, (17) 



which happens to be identical to the Rains bound (14) in 
the whole region AYCD. So the upper and lower bound 
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FIG. 4: A close-up of the non-additive OO-invariant states in 
Figure 2, for the purpose of calculating the minimal convex 
extension to the AREEP. In region AYCD, the minimal con- 
vex extension only depends, afnnely, on / (eq. 17). In region 
BCY, the minimal convex extension is afBne along the lines 
depicted here (given by eq. 18). 



end point (segment XY). By convexity of Er and of the 
Rains bound, and by the fact that the Rains bound is an 
upper bound on Er and the tangent a lower bound, it 
follows that both Er and the Rains bound must coincide 
with this tangent and, hence, be affine along each of the 
lines (18). We conclude that Er is equal to the Rains 
bound also in the remaining region YCB. 

F. Summary of results 

We finalise the calculation of Er on the OO-invariant 
states by summarising all the results obtained for the 
different regions in the following table. Figure 5 shows a 
contour plot of Er for the case d = 3. 

Region Er 
PPT 

A E R , (9), with s = 1+(d ~y/~* and s = 1 

B \ AYCD E R , (9), with s = and s = j^j 
C \ CYB E R , (9), with s = and s = 1 
AYCD eq. 14 

CYB affine along lines (18) between YX and BC 
YX eq. 15 

Furthermore, the Rains bound is equal to Er in any of 
these regions. 

VII. CONCLUSION 



equal each other within this region and, hence, Er is 
equal to the Rains bound in AYCD. 

The situation for the remaining triangle YCB is some- 
what more complicated. To calculate Er we consider a 
set of straight lines connecting points on the line segment 
BC with points on the segment XY and given by 

- p(d 2 -2 + (rf-2)p) 

f = ~Pf+ 2 + d(p-2)-2p ■ (18) 

These lines are parameterized by p, which runs from 

to =g. Recall that the line XY is given by / = (l + /)d/2 

and BC by / = 3 - 4/d + (d - 1)/. 

On the line segment XY, the Rains bound is given 
by (15). On the segment BC, and in fact to the right 
of it as well, the Rains bound is equal to Er = Er 
and is given by (9) with s — and s = 1 (region C of 
Figure 1). Moreover, this formula holds for all points 
on the lines (18) within the additivity region, allowing 
for the calculation of the derivative of the Rains bound 
along the lines (18). Doing this in the points on the 
additivity border BC yields the result that, for every line 
(18), the tangent to the Rains bound at the start point 
(on segment BC) touches the Rains bound again at the 



In this paper, we have considered the calculation of the 
AREEP Er for the class of OO-invariant states, gener- 
alizing the results of [10], which dealt only with the class 
of Werner states. This has been achieved using four ba- 
sic ingredients: properties of the REEP Er, properties of 
the Rains bound R (1), and a deep connection between 
these two quantities Er and R. The final cornerstone 
of the calculation is the symmetry inherent in the OO- 
invariant states [6]. 

The relevant properties of the REEP are that it is an 
additive entanglement measure in a large region of state 
space [4] and that the AREEP is convex everywhere [9] . 
This convexity allows us to use the "minimal convex ex- 
tension" construction as a lower bound. 

We have shown here that the Rains bound is also con- 
vex and continuous, and that the calculation of it can be 
reformulated as a convex optimization problem, which 
implies, by the way, that this problem can be solved effi- 
ciently and does not suffer from multiple local optima. 

We have also made explicit the techniques that were 
already employed in [10] implicitly, resulting in Lemma 
4. This Lemma shows that there is a deep connection 
between the AREEP and the Rains bound and seems to 
suggest that both regularise to the same quantity [7] . Un- 
fortunately, in its current form, the Lemma is weakened 
by the additional requirement on the states er, over which 
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FIG. 5: Contour plot of the AREEP E'g for the OO-invariant 
states, parameterized by / and / (d = 3). Superimposed on 
this plot are the lines separating the different regions defined 
in the text (regions A, B and C, the PPT set, the set of 
additive states and the regions AYCD and CYB). 



the Rains bound is minimized, that the quantity \<r T2 



should be positive. We have coined the term binegative 
states for those states that violate this requirement and 
we have made some initial investigations into the ques- 
tion of their existence. Specifically, we showed that for 
the case of OO-invariant states, a is not binegative, so 
that the Lemma can be used here at full strength. If it 
turned out that the extra requirement can always be re- 
moved, in one way or another, then the Lemma could di- 
rectly be used to prove Rains' suggestion that E^ = R°° . 

For the time being, we have been able to show that 
at least for ^-additive states p, the Rains bound and 
the REEP are equal (and, of course, also equal to their 
regularised versions). 

Using these results, we have calculated the AREEP 
for OO-invariant states and it followed as a by-product 
of the calculation that the Rains bound is identical to 
Eft for the OO-invariant states. 

This last result could be taken as a hint that the Rains 
bound might be additive everywhere, in contrast to Er. 
If this were true, then this would imply that the AREEP 
is precisely equal to the (non-regularised) Rains bound 
and, furthermore, that it can be calculated efficiently. 

This work has been supported by project GOA- 
Mefisto-666 (Belgium), the EPSRC (UK), the European 
Union project EQUIP, the Deutsche Forschungsgemein- 
schaft (Germany) and the European Science Foundation. 
KA wishes to thank J. Eisert, M. Plenio, F. Verstraete, 
J. Dehaene and E. Rains for fruitful discussions and re- 
marks. 



[1] P. Horodecki and R. Horodecki, Quant. Inf. Comp. 1, 45 
(2001); R.F. Werner, Quantum information - an intro- 
duction to basic theoretical concepts and experiments, in 
Springer Tracts in Modern Physics 173 (Springer, Hei- 
delberg, 2001); W.K. Wootters, Quant. Inf. Comp. 1, 27 
(2001); M.A. Nielsen and G. Vidal, Quant. Inf. Comp. 1, 
76 (2001); M.B. Plenio and V. Vedral, Cont. Phys. 39, 
431 (1998). 

[2] R.F. Werner, Phys. Rev. A 40, 4277 (1989). 

[3] V. Vedral, M.B. Plenio, M.A. Rippin, and P.L. Knight, 

Phys. Rev. Lett. 78, 2275 (1997); V. Vedral and M.B. 

Plenio, Phys. Rev. A 57, 1619 (1998); V. Vedral, M.B. 

Plenio, K.A. Jacobs, and P.L. Knight, Phys. Rev. A 56, 

4452(1997). 

[4] E. Rains, Phys. Rev. A 60, 179 (1999). Erratum: Phys. 



Rev. A 63, 019902 (2001). 
[5] E. Rains, IEEE Trans. Inf. Theory 47(7), 2921-2933, 
2001. 

[6] K.G. H. Vollbrecht and R.F. Werner, Phys. Rev. A 64, 

062307 (2001). 
[7] E. Rains, private communication (2001). 
[8] M. Ohya and D. Petz, Quantrum Entropy and its Use, 

Springer- Verlag, 1993. 
[9] M. Donald, M. Horodecki and O. Rudolph, Lanl e-print 

0105017. 

[10] K. Audenaert, J. Eisert, E. Jane, M.B. Plenio, S. Virmani 
and B. De Moor, Phys.Rev.Lett. 87, 217902, 2001. 

[11] M. Donald and M. Horodecki, Phys.Lett. A 264, 257-260 
(1999). 



